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We argue that the observed antiproton production in heavy-ion coUisions at CERN-SpS energies 
can be understood if (contrary to most sequential scattering approaches) the backward direction 
in the process pp nvr (with n=5-6) is consistently accounted for within a thermal framework. 
Employing the standard picture of subsequent chemical and thermal freezeout, which induces an 
over-saturation of pion number with associated chemical potentials of /^^ ~ 60-80 MeV, enhances the 
backward reaction substantially. The resulting rates and corresponding cross sections turn out to be 
large enough to maintain the abundance of antiprotons at chemical freezeout until the decoupling 
temperature, in accord with the measured p/p ratio in Pb(158AGeV)-|-Pb collisions. 
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Over the last decade remarkable progress has been 
made in the understanding of the dynamics of strong in- 
teractions probed through (ultra-) relativistic heavy-ion 
collisions. Although the main challenge of an unambigu- 
ous identification of the QCD phase transition to the 
Quark-Gluon Plasma (QGP) persists, we have greatly 
advanced our knowledge on properties of highly excited 
hadronic matter close to the expected phase boundary. 
A variety of collective phenomena has been observed in- 
dicating that the produced systems have indeed reached 
macroscopically large sizes, justifying the use of equilib- 
rium techniques such as thermo- and hydrodynamics. 

One of the important results that will be used be- 
low is that the final-state hadron abundances, including 
antibaryons, can be rather accurately characterized by 
the so-called chemical freezeout stage Q with a com- 
mon temperature Tch and baryon chemical potential 
/i^, the specific values depending on collision energy 
(note that a precise description of all hadron species re- 
quires corrections to an ideal gas ensemble, e.g., excluded 
(eigen-) volumes ||^,^ to mimic short-range repulsions, or 
'strangeness suppression' factors ||^ . Such corrections are 
not important for our subsequent analysis and will be ne- 
glected. For a contrasting view of hadron production in 
heavy-ion collisions, see, e.g., ref. (s)). 

At SpS energies, the chemical freezeout is clearly dis- 
tinct from the thermal one (with an associated temper- 
ature Tth < Tch), from where on the particles stream 
freely to the detectors. This follows from the kinetics of 
hadronic reactions i.e., a significant difference be- 

tween elastic and inelastic collision rates at low relative 
energies, and has been confirmed by several experimen- 
tal evidences (see, e.g., ref. In central Pb-|-Pb colli- 
sions, e.g., a nucleon at midrapidity is elastically rescat- 
tered on average about 10-15 times, but less than once 
inelastically 0. Also, as shown in most of the collec- 
tive flow effects at SpS are generated in between the two 
freezeouts, and their observation leaves no doubt about 
the existence of such an intermediate stage. 

Another consequence is that abundances of secondary 
mesons (pions, kaons, etc.) are not subject to signifi- 
cant changes when the system evolves from Tch to Tth. 
Again, this is conceivable as the rates for number chang- 



ing reactions are too small to maintain chemical equilib- 
rium on the time scales of the hadronic fireball lifetime, 
t ~ 10 fm/c. Since, on the other hand, elastic scat- 
tering {e.g., TTTT —tp—i- tttt) is still effective, thermal 
equilibrium is approximately maintained. In a statisti- 
cal mechanics language such a scenario entails additional 
conservation laws, which can be implemented via effec- 
tive (pion-, kaon-, etc.) chemical potentials to guarantee 
fixed particle numbers. For pions (kaons) typical values 
of = 60-80 MeV (/^^ = 100-130 MeV) are reached 
at Tth = 110-120 MeV, see @,|,|. The ensuing 'over- 
saturation' of the pion phase space is to play the key role 
in what follows. 

The situation for antibaryons is different from mesonic 
secondaries as the pertinent inelastic (or annihilation) 
cross section is not small. At the relevant (thermal) 
energies in the comoving frame of collective expansion, 
y/s — 2{mN + E*j^) ~ 2.3 GeV, one has crpp_>„^ ~ 50 mb. 
Taking an average baryon density of qb = 0.75go {go — 
0.16 fm"'^) around T — 150 MeV in the course of the 
hadronic evolution and a typical (anti-) proton velocity of 
vth = p/Etot = 0.56c (from E^'' = {3/2)kT ~ 225 MeV), 
we obtain for the chemical equilibration time scale 
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Tch 



C^ann QB Vth 



3 fm/c . 



(1) 



This is well below the fireball lifetime in the pure 
hadronic phase of Thad — 7 fm/c cf. Fig. [l| in 

other words, only a fraction of Gxp\—- Thad /Tch] — 10 % 
does not rescatter towards thermal freezeout. Thus, in 
spite of the good agreement of the final p yields with 
the standard chemical freezeout predictions, antibaryon 
production ought to be reconsidered. 

Naively one might expect most of the antiprotons to 
be annihilated, and various transport calculations {e.g., 
ARC § and UrQMD ^) have indeed been unable to 
account for the measured number, falling short by signifi- 
cant factors. Consequently, speculations have been raised 
that the puzzle could be resolved by certain in-medium 
effects leading to either a reduction of the annihilation 
rate p|] or an enhanced production pOl . 
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1. Chemical relaxation time (full line) for antipro- 



tons in central Pb + Pb collisions at SpS energies employ- 
ing an annihilation cross section of a ann~ (4:0-50) mb and a 
baryon-density profile as obtained within an isentropically ex- 
panding fireball model ^ . The dashed-dotted curve indicates 
the remaining hadronic fireball lifetime till thermal freezeout. 

A generic problem with transport/cascade simulations 
is that, although multi-body resonance decays are in- 
cluded, the inverse reactions are not. As is well-known, 
this violates detailed balance and, in principle, prevents 
the simulations from reaching the proper thermodynamic 
limit |pl| . The discussions are typically focused on three- 
body reactions; in the present context it is even 5- or 
6-pion collisions which are most relevant for producing 
baryon-antibaryon pairs. A widespread belief is that 
those reactions have insignificant rates under the con- 
ditions probed in the hadronic stages of heavy-ion colli- 
sions. However, as we will argue below, this is not the 
case; inverse reactions have to be addressed and are, in 
fact, capable of explaining the observed antiproton yield 
(although the p enhancement in nucleus- nucleus over p-p 
collisions thus looses its proposed direct relation to QGP 
formation, the latter is by no means excluded by the sub- 
sequent arguments). 

The expression for the thermal reaction rate for the 
process pp ^ mr can be written as (see also Ref. ||l^ ) 

Uth = Jd% d% (f'K, ■ --d^'K^ S'^^\Ktot) \Mn\^ 

x{zp Zp exp 



'^]-z^ exp[-2^— ] } , (2) 
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where A1„ denotes the invariant scattering matrix ele- 
ment (which, of course, is identical for the back- and for- 
ward direction due to time-reversal invariance of strong 
interactions), d'^k^ = d^kx/i^n)^ the phase space inte- 



grations and Kt, 



kn k-TT 



the total 



four-momentum; Zx — e^^/-^ are the fugacities of par- 
ticle species x (in Boltzmann approximation), and n is 



the number of pions produced in a pp annihilation at 
a given energy (to be discussed in more detail below). 
For a nonvanishing net nucleon number and in chemical 
equilibrium, one has > 0, /i^ = —^n and /i^r — 0. 
Under SpS conditions, typical values at chemical freeze- 
out are § (r^/^^i^) ~ (170,260) MeV, which results in 
an antiproton-to-proton ratio of 
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(Ep + tiM)/T] 



exp[-2/^jv/T] = 4.7% 



(3) 



consistent with the experimentally measured value of 
(5.5 ± 1) % [^3[[l|]. Note also that at chemical freezout 
ZpZp = 1, = 1, and the forward and backward rates in 
eq. (H) are simply equal (the Boltzmann factors in both 
terms contain the same total energy). This implies, e.g., 
that cascade simulations ignoring the back-reaction at 
comparable pion densities [qt,: ~ 0.2 — 0.25 fm"-^) cannot 
give a proper account of the antiproton production. 

The evolution of the system towards thermal freezeout 
can be constructed using a thermal fireball model which, 
based on isentropic expansion, leads to (Tj^i,^^) ~ 
(120, 415) MeV. At first sight, this gives a p/p ratio 
of exp[— 2^Ar/r] ~ 0.1 %, a factor of ~50 below the 
experimental result. However, this estimate is lacking 
an important ingredient. A correct statistical treatment 
including pion-number conservation forces us to intro- 
duce individual chemical potentials/fugacities for anti- 
/protons and pions. Insisting on equilibrium for the re- 
action in question, the following relation holds: 
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(4) 



Since antibaryons constitute only a small fraction of the 
produced secondaries, they have an insignificant impact 
on the evolution of nucleon- and pion-chemical potentials. 
Therefore, eq. (Q) provides an estimate of the antiproton 
fugacity at thermal freezeout. The crucial point here is 
that pion over-saturation will generate a strong enhance- 
ment of the back-reaction n7r — > pp. 

For a more quantitative assessment it is important to 
have a rather accurate determination of the pion multi- 
plicity distributions in pp annihilation. A nice discussion 
of its systematics has been given in ref. |l|l, which we will 
follow here to incorporate the empirical knowledge. For 
pp annihilation at rest (^/s — 2mjv), the data can be rep- 
resented by a least-square fit to a Gaussian probability 
distribution (normalized to one). 
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with a mean (n) = 5.02 and width cr — 0.90 (see also 
ref. |[T^). At larger cm. -energies both the average num- 
ber and width increase. For the former, linear |17| and 
logarithmic |lq| dependencies have been proposed. 



(n) = ci + C2 .s°-5 
(n) = Ci + C2 ln(s) 



(6) 
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(s in [GeV^]), which both reproduce the measured multi- 
phcities up to at least -y/s = 5 GeV using the parameters 
ci = 2.6 ± 0.5, C2 = (1.3 ± 0.2) GeV-^ and ci = 2.65, 
C2 — 1.78, respectively. For the latter, the energy depen- 
dence of the width has also been given as [E8| 



= 0.174 



„0.2 



(7) 



For our application in a thermal environment at T = 
150 MeV (implying ^ = 2.33 GeV) we fix (n) = 5.65 
(in accordance with eq. (^)) together with a 10% increase 
in a (as suggested by eq. (0)) to extract discrete weights 
Wn = P{n; {n),a). The averaged pion-fugacity enhance- 
ment factor then follows as 



E 

n=2 



exp[n/x^/r] 



(8) 



where Umax = 9 for any practical purpose. Insert- 
ing now thermal freezeout values Tth — 120 MeV and 
/ij/* ~ 65 MeV (as arising in a thermal fireball model [^), 
yields (z") — 25. This entails a large enhancement of the 
antiproton-to-proton ratio, from 0.1% to 2.5%. In fact, 
owing to the high power of the pion fugacity, slightly 
larger chemical potentials of = 75-80 MeV result in 
an enhancement factor of 42-54, rendering the pertinent 
p/p-ratio in line with the observed (chemical freezeout) 
value, cf. Fig. |2[ 
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FIG. 2. Antiproton-to-proton ratio as a function of (de- 
creasing) temperature in an isentropically expanding fireball. 
The dashed curve represents the naive ratio, exp[— 2^jv/r], 
whereas the full curves are for finite pion chemical potentials 
indicating uncertainties as discussed in the text. The experi- 
menatl data point is from Ref. [ p^ . 

Such slightly increased values for the pion chemical 
potential close to thermal freezeout can indeed be easil 
argued for. Within the thermal fireball model of ref. 
elastic ttN B scattering {B: baryonic resonances up 
to niB — 1.7 GeV) was assumed to be frequent enough to 



maintain (relative) chemical equilibrium for the occupa- 
tion of the excited baryonic states. However, with typical 
corresponding cross sections of cTt^n^b — 15-30 mb [ p^ , 
this might not be fully justified anymore for the last few 
fm/c prior to thermal freezeout. Consequently, a larger 
fraction of the pion number resides in explicit pionic de- 
grees of freedom rather than in excited resonances, which 
translates into an effectively larger ^lj^. 

Let us finally comment on implications of our obser- 
vations for RHIC. Close to the expected chemical freeze- 
out the pion density is very similar to SpS conditions. 
Thus the rate of producing antiprotons through multi- 
pion annihilation per unit time and volume is essentially 
the same in both cases. The crucial difference is, how- 
ever, that the total density of antiprotons is much larger 
around midrapidity at RHIC due to substantially smaller 
baryon chemical potentials. More quantitatively, using 
typical thermal model estimates with p^* ~ 0.2 go 
shortly after chemical freezeout (further reduced there- 
after), one obtains t^^^ ~ 11 fm/c. With the life- 
time of the hadronic phase at RHIC being compara- 
ble to that at SpS energies, chemical equilibrium in the 
pp <-> nTT reaction cannot be maintained until thermal 
freezeout (also, the emerging pion oversaturation is less 
pronounced in a baryon-poor regime). The observed 
antiprotons at RHIC should therefore mostly originate 
from earlier stages, corresponding to the standard hadro- 
chemical freezeout in the vicinity of the phase boundary. 
Nevertheless, our time scale estimate indicates that even 
under RHIC conditions, antibaryon annihilation will be 
partially compensated by the inverse reactions. 

To summarize, we have analyzed the p/p-ratio at SpS 
energies employing a thermal approach. So far this ob- 
servable has been difficult to understand within, e.g., 
transport models which only included the annihilation 
channel, causing doubts whether the latter is actually 
active, or unconventional mechanisms for enhanced pro- 
duction need to be invoked. We have shown, however, 
that the 'puzzle' can be resolved in a rather standard 
statistical-mechanics framework upon inclusion of the in- 
verse process of multipion scattering into pp pairs, which 
can be supported until thermal freezeout. Our main 
ingredient was that effective pion-number conservation 
generates pion over-saturation at the later stages of a 
heavy-ion collision, as described by the build-up of appre- 
ciable pion chemical potentials. Raised to a large power 
(n ~ 6) the corresponding pion fugacities sustain a high 
antibaryon fraction, thus counter-balancing the loss from 
BB annihilation. This mechanism also complies with the 
measured centrality dependence being essentially con- 
stant, as to be expected from a hadro-chemistry varying 
little with impact parameter (for sufficiently peripheral 
collisions the applicability of thermal model analyses, of 
course, ceases and p production, normalized to the num- 
ber of participant nucleons, approaches its value in p-p 
collisions, which lies about 30% below the one in central 
nucleus- nucleus reactions Q). 

Finally we should note again that our findings are not 
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contradictory to an earlier QGP formation. They rather 
give further support to the equiUbrium concept of subse- 
quent chemical and thermal freezeout stages, which has 
already proven successful in the explanation of a vari- 
ety of observables, such as hadron abundances, collective 
flow, HBT radii and electromagnetic radiation. 

This work was supported by the U.S. Department of 
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[1] P. Braun-Munzinger, J. Stachel, J. P. Wessels and N. Xu, 
Phys. Lett. B365, 1 (1996); 

P. Braun-Munzinger, I. Heppe and J. Stachel, Phys. Lett. 
B465, 15 (1999). 
[2] G.D. Yen and M.I. Gorenstein, Phys. Rev. C59, 2788 
(1999). 

[3] J. Rafelski and J. Letessier, ^ep-ph/9903018 
[4] E. Shuryak, Phys. Lett. B207, 345 (1988); Pliys. Rev. 
D42, 1764 (1990); 

P. Gerber, H. Leutwyler and J. L. Goity, Phys. Lett. 
B246, 513 (1990); 

H. Bebie, P. Gerber, J.L. Goity and H. Leutwyler, Nucl. 
Phys. B378, 95 (1992). 
[5] CM. Hung and E. Shuryak, Phys. Rev. C57, 1891 
(1998). 

[6] R. Stock, Nucl. Phys. A661, 282c (1999). 
[7] H. Sorge, Phys. Rev. C52, 3291 (1995). 
[8] R. Rapp and J. Wambach, Eur. Phys. J. A6, 415 (1999). 
[9] Y. Pang, D.E. Kahana, S.H. Kahana and H. Crawford, 
Phys. Rev. Lett. 78, 3418 (1997). 
[10] M. Bleicher, M. Belkacem, S.A. Bass, S. SofT and H. 

Stocker, Phys. Lett. B485, 133 (2000). 
[11] M. Belkacem et ai, Phys. Rev. C58, 1727 (1998). 
[12] P. Koch, B. Miiller and J. Rafelski, Phys. Rep. 142, 167 
(1986). 

[13] M. Kaneta for the NA44 collaboration, Nucl. Phys. 

A638, 419c (1998). 
[14] F. Sikler for the NA49 collaboration, Nucl. Phys. A661, 

45c (1999). 

[15] C.B. Dover, T. Gutsche, M. Maruyama and A. Faessler, 

Prog. Part. Nucl. Phys. 29, 87 (1992). 
[16] W. Bliimel and U. Heinz, Z. Phys. C67, 281 (1995). 
[17] P. Gregory et al, Nucl. Phys. B102, 189 (1976). 
[18] R. Stenbacka et al, Nuovo Cim. 51A, 63 (1979). 
[19] Particle Data Group, D.E. Groom et al, Eur. Phys. J. 

C15, 1 (2000). 

[20] S. Bass et ai, Proc. of Quark Matter '99, Nucl. Phys. 
A661, 205c (1999). 



4 



